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ABSTRACT 

Pluto and its five known satellites form a complex dynamic system. Here we explore where 
additional satellites could exist exterior to Charon (the innermost moon) but interior of Hydra 
(the outermost). We also provide dynamical constraints for the masses of the known satellites. 
We show that there are significant stable regions interior of Styx and between Nix and Kerberos. 
In addition, we show that coorbitals of the known small satellites are stable, even at high 
inclinations, and discuss mass constraints on undiscovered satellites in such orbits. 


Subject headings: Kuiper belt objects: individual(Pluto) — planets and satellites: individual(Pluto) — 
planets and satellites: dynamical evolution and stability — celestial mechanics 


1. Introduction 

The discovery of Pluto’s major satellite Charon 
by IChristv &: Harrington ( 1978 1 made Pluto- 
Charo n the first known n ear-equal mass ratio 
(8.2:1, Brozovic et al. 2015ll in the solar system, 
and Charon remains the largest known satellite 
in the Kuiper Belt. Pluto and Charon orbit their 
common barycenter in near-circular orbits with 
a period of 6. 3872 days and an ave r age distance 
of 19,596 k m ( Brozovic et al. 201, 'ill . Stern et al. 
(|l99lLll994l l searched for additional satellites from 


the ground and with the Hubble Space Telescope 
(HST), and while they did not find any satel¬ 
lites, the later study did discover the Charon In¬ 
stability Strip, which destabilizes any any small 
prograde satellites withi n the 2:1 mean mot ion 
resonance with Charon. Weaver et al. ( 2006ll fi¬ 


nally discovered two small satellites, Nix and Hy- 


2011 (Showalter et al. 

2011 ) 

Showalter et al. 20121). 

The 


four small satellites are closely packed together 
and are stro ngly perturbed by b oth each other 
and Charon ( Youdin et al.l 12012 ). In addition, 


Giuliatti Winter et al. ( 201(1 201,31 2014 ) have 
shown that small satellites can be stable between 
Pluto and Charon, even at high eccetricities and 
inclinations. 

NASA’s New Horizons spacecraft will fly 
through the Pluto system in July 2015, passing 
within 12,500 km of Pluto and 1,800 km of the 


Pluto-Charon L 3 point ( Guo fc Farauhai 20081) . 
Because of this close passage, the spacecraft will 
perform an intensive search for new satellites and 
rings of Pluto to make sure that there are no 
objects currently producing dust particles which 
could damage the spacecraft. Though the like¬ 
lihood of a new satellite being hazardous to the 
spacecraft can be shown to be very low, the search 
for new satellites on approach will be the most 
complete survey of small Pluto satellites possi¬ 
ble, far exceeding what can currently be done 
with ground-based and Earth-orbit observatories. 
Here, we show where it is dynamically possible to 
insert new satellites with a small non-zero mass 
in the Pluto system and speculate which of is¬ 
lands of stability are likely to be populated with 
currently-unseen satellites or rings. 
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2. Numerical Methods 

To simulate the Pluto system, we implemented 
an n-body integrator in C++ with a Python inter¬ 
face. This integrator is based on the 12(10) Runge- 
Kutta-Nystrom meth od in ACM Algorithm 670 


([Brankin et al.l 1 19891) , and our implementation 
is available onlin^il- This method is optimized 
for systems with a few strongly-interacting ob¬ 
jects and no non-gravitational nor tidal dissipa¬ 
tion forces. We verified this method against an 
independent high-precision Bulirsch-Stoer inte¬ 
gration of the Pluto system for 1000 years, and 
found agreement between the two integrations to 
within 10 meters (D. P. Hamilton 2014, private 
communication). For this implementation, we set 
a maximum distance from the barycenter of 10® 
km, i.e., the simulation would be stopped if any 
object passed beyond this distance. In addition, 
we assigned effective collision radii for the known 
and new small satellites of 500 km, so that the sim¬ 
ulation would fail if any two small bodies passed 
within 1000 km of each other. While in reality 
such an encounter is not close enough to be an 
impact, we believe it is close enough to call into 
question the long-term stability of those two ob¬ 
jects (the Hill radii of the known small satellites 
are all <20 km). 

2.1. Known Satellite Uncertainties 


Brozovic et al.l (12015^ is the best published esti¬ 


mate of the masses and orbits of the known Pluto 
satellites. However, while they provide mass un¬ 
certainties, the mass and state vectors of the satel¬ 
lites are coupled and certain masses are only pos¬ 
sible with certain states. We explored this param¬ 
eter space by creating a dense discretely-sampled 
probability distribution function (PDF). We con¬ 
structed an approximate log-probability function 
for a given state using the measure ment uncer¬ 


tainty of the HST astrometric poin t s iniBuie et al 
( 20061 201, “ih and Brozovic et al. ( 2015 ). enforc¬ 


ing a minimum lifetime of 10® Charon periods, 
and rejecting solutions where any of the ob¬ 
jects passed within 1000 km of each other or 
one of the known satellite’s barycentric eccen¬ 
tricity exceeded 0.1. Since the uncertainty in 
Charon’s orbit is much smaller, we used the best- 


ht Pluto and Charo n masses and states from 
Brozovic et al. ( 20151) for all the simulations. We 
then fed this function into the “emcee” Affine 
Invariant Markov chain Monte Carlo (MCMC) 


Ensemble sampler (jForeman-Mackev et al.l 120131) 


to build the PDF. The result was an ensemble 
of probability-weighted possible known-satellite 
mass/state systems from which we could draw. 
Because we enforced a minimum lifetime, our me¬ 
dian masses fo r the known sa t ellites were slightly 
different than Brozovic et al.l ( 2015h . but within 
b oth their error bars a nd the stability constraints 


of lYoiidin et al.l (|2012l) . 


^https://github.com/ascendingnode/PyNBody 


2.2. Simulations of New Satellites 

We simulated unknown/unseen objects in the 
Pluto system adding an object with a specified 
non-zero mass and randomized state vector less 
than a given inclination, verifying that it does not 
make the residuals of the known satellites signif¬ 
icantly worse (<5% increase in y^), and making 
sure that it has a lifetime greater than 10® Charon 
orbits. We assumed the new satellites had a grav¬ 
itational mass of 10“® km® s“^, about 10% of the 
mass of our best fit for Styx. To create a random 
state vector, we first generated two random unit 
vectors for bar ycentri c posit ion and velocity using 
the method in iKno 3 (Il970l) and made sure that 
their cross product (i.e. the direction of the satel¬ 
lite’s angular momentum vector) was angled less 
than a given inclination from the Pluto-Charon 
orbit pole. For the final simulations, we gener¬ 
ated cases with <10° inclination (to fully explore 
coplanar satellites) and a general case with un¬ 
restricted inclination. We then generated a uni¬ 
formly random magnitude for the position vec¬ 
tor between 20,000 km (just outside the orbit of 
Charon) and 100,000 km (beyond the orbit of Hy¬ 
dra) . We calculated the circular orbital velocity at 
this distance from the barycenter and then gener¬ 
ated a random velocity between 0 and 1.5 times 
that circular velocity. We also made sure that 
the period of this orbit was less than seven times 
that of Charon, as orbits beyond that point are 
not significantly influenced by the known small 
satellites. For each state vector that passed these 
tests, we next added to the known objects of the 
Pluto system and checked that the resulting 
w as smaller than 1 .05 tim es the y^ of the solution 
in iBrozovic et al] ( 2015 ). Finally, objects which 
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survived for the ten-year span of observations were 
integrated forward for 10® Charon orbits to verify 
their medium-term stability against close encoun¬ 
ters with the known satellites. Any test objects 
that passed within 1000 km of a known satellite 
were deemed unstable and rejected. 

3. Results 

Our assumed mass distribution for the known 
satellites is shown in Figure [TJ and compared to 
previous estimates in Table [TJ Figure [2] shows the 
stable near-coplanar barycentric orbits, while Fig¬ 
ure O shows the orbits that were unrestricted in 
inclination. In Figures [5] and [31 the blue points 
are all the orbits, while the red points are the 10 % 
with the lowest astrometric residuals. 


3.1. Refined Known Satellite Masses and 
Implications 


We were able to obtain much tighter con¬ 
straints on the masses of the known satellites 
than previous studies by enforcing stability for 
10® Ch aron orbits (1747 year s), rather than 200 


years in iBrozovic et al.l (1201511. O ur largest diver¬ 


gence from Brozovic et al. ( 2015h was that by re¬ 
quiring Styx to be both stable and have a positive, 
non-zero mass, we obtained a smaller mass for Nix 
(Nix is the largest perturber of Styx after Charon). 
The smaller mass of Nix then required the mass 
of Hydra to grow, in order to place the same per¬ 
turbations on Kerberos. Our Kerbe ros m ass is 


almost the same as Brozovic et al 


and 


suppo rts the prediction in Ishowalter fc Hamilton 


(12015 1 that Kerberos 1V=25.95. IShowalter et al 
201111 must have a significantly darker albedo 
than Nix and Hy dra (V=23.38 and V=22.93, 
Weaver et al.M2006ll . The apparently dark surface 
of Kerberos is curious given that the small satel¬ 
lites can exchange surface material by sw eeping 
up each other’s low-velocity im pact ejecta (ISternI 
2008 : Porter fc Grundv 201 Sll . Our Hydra/Nix 
mass ratio (Ril.5:l) is different to the mass ratio 
obtained by using the flux ratio and otherwise 
assuming equal properties (~ 1 . 8 : 1 ), but they are 
close enough that shape e ffects may be enough to 
accou nt for the difference ([Showalter fc Hamilton 
201511 . 


3.2. Possible Orbits of Near-Coplanar 
Satellites 


As shown in Figure HJ we found a great diver¬ 
sity of potential orbits for low-inclination unseen 
satellites which were both stable to the known 
satellites and did not adversely affect observed 
astrometric residuals. The most stable region is 
just interior of Styx, centered just interior of the 
3:1 mean motion resonance (MMR) with Charon. 
Similar, slightly less stable regions are between 
Styx and Nix (around the 11:3 MMR) and Nix 
and Ke rberos (aroun d the 9 :2 MMR). As identi¬ 
fied bv lYoudin et al.l ( 2012 1. the region between 
Kerberos and Hydra is very unstable. Orbits be¬ 
yond the 13:2 MMR were effectively unconstrained 
by the known satellites, while all prograde orbits 
in terior of the 2:1 M MR were unstable, as shown 
in IStern et al.l ( 1994ll . 


Regardless of period, we found that the num¬ 
ber of stable orbits increases with inclination at 
a rate faster than cos(inchnation). This is due 
to lower-inclination orbits interacting more often 
and strongly with the known objects (which are all 
within 1 ° inclination), thus producing a more no¬ 
ticeable effect on the residuals. Orbits interior of 
Styx are stable at the lowest inclinations without 
producing a noticeable effect on the residuals. 


We found that coorbitals of all the known satel¬ 
lites are also stable. In the rotating frame of the 
known satellite they are coorbiting with, we found 
that most of these sate llites are in low-ecce n tricity 
“horseshoe orbits” ( Murray fc Permed 119991 ) 
with large libration amplitudes, though some 
(mainly Hydra coorbitals) were in “tadpole or¬ 
bits” librating around the L 4 or L 5 Lagrange 
points. Nix, Kerberos, and Hydra coorbitals with 
our assumed mass of 10 “® km® s“^ had a negative 
effect on the residuals, so any real coorbitals must 
be significantly lower mass than this. However, 
Styx coorbitals with our assumed unknown satel¬ 
lite mass (ssl0% Styx) had a much smaller effect 
on the residuals, though this may be because we 
had much fewer astrometric points for Styx than 
the other satellites, so the problem is less well 
constrained. 
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3.3. Possible Orbits of High Inclination 
and Retrograde Satellites 

As seen in Figure |31 we found stable high in¬ 
clination orbits. At high (near-polar) inclinations, 
the region between Nix and Kerberos is most sta¬ 
ble, followed by the region between Kerberos and 
Hydra. Orbits interior of Nix at high inclinations 
are generally unstable. 

There is also a huge stability plateau for retro¬ 
grade satellites around the -7:3 MMR. These or¬ 
bits are stable to Charon and far enough interior 
of Styx that they do not affect the residuals. In 
addition, there are curious Kerberos and Hydra 
coorbital structures at ?»140° inclination. 


4. Discussion 


We have shown that there are stable niches 
between Charon and Hydra at low inclination. 
If there are any presently-undiscovered low- 
inclination Pluto satellites between Charon and 
Hydra, it is likely they formed in the same 
way as the known satellites, which are gener¬ 
ally assumed to be remainin g fragments from 
the Charon-forming imp act ( Stern et"^ l2006t 
iKenvon fc Bromic 2014 and citations therein). 
The most likely region for an unseen satellite is 
therefore in a low-inclination, low-eccentricity or¬ 
bit with a period between 16 and 19 days. Satel¬ 
lites up to 10% of the mass of Styx can stably 
exist and not adversely affect the orbital residu¬ 
als of the known satellites here. Low-inclination, 
low-eccentricity satellites with periods of around 
24 and 29 days are also possible, but would need 
to be very low mass (<1% Styx) to avoid affecting 
the observational residuals. 

This is in line with the detection limits from 
previous HST-based Pluto satellite searches, 
which placed strict limits on the brightness of 
objects between Nix and Hydra, but were much 
less se nsitive interior of Sty x due to the glare from 
Pluto ( Steffi fc Sternl[200^ . Since we were able to 
determine a dynamical mass for Styx of 10“"^ km^ 
s“^, our test mass of 10“® km^ s“^ is just small 
enough to have plausibly passed detection previ¬ 
ously. Our test particles could have a noticeable 
perturbing effect on the real satellites for any or¬ 
bit between Styx and Hydra, suggesting that the 
dynamical limits to the size of any low-inclination 
satellites in this region are even more restrictive 


than the limits from HST-based searches. 

NASA’s New Horizons spacecraft will perform 
a deep search for new satellites on approach to 
Pluto, using the Long-Range Reconnaissance Im¬ 
ager ( LORRI) with a sen sitivity of greater than 


V=17 (ICheng et al.ll2n08h in order to assess the 


hazard to spacecraft from ejecta dust produced 
currently-undiscovered satellites. The final hazard 
observations take place at a distance of less than 
0.1 AU from the Pluto barycenter; at this distance, 
LORRI is sensitive to objects down to Hy > 14, 
or a few kilometers across. New Horizons could 
therefore discover satellites much smaller than our 
mass constraints. 

We also found several classes of stable high- 
inclination and retrograde satellites. Some of 
these orbits, particularly those between Nix and 
Kerberos could exist at both low and high incli¬ 
nations, meaning a former low inclination satellite 
in this period range could be perturbed by an im¬ 
pact to higher inclination. However, such a lucky 
impact is highly unlikely. 

The large stable retrograde region interior of 
Styx could not easily by populated by remnants 
from the Charon-forming impact, but could in¬ 
stead ho st captured satellite s . Binary-exchange 
captu re ( Funato et al.l 120041 : lAgnor fc Hamilton 
20061 ) would happen when a binary KBO has 
very close encounter with the Pluto system, and 
one part of the binary is captured, while the 
other part is ejected from the system. Such 
an event is very unlikely today, but would have 
been much more likely in the giant planet migra¬ 
tion era, when Triton is presumed to have been 
captured around Neptune in exactly this fashion 
( Agnor fc HamiltonI l2006^ . Such a lucky stable 
capture is very unlikely, but if it happened, it 
would almost certainly be to a retrograde orbit 
with a period between 13 and 19 days. 

We did not perform any direct simulations of 
rings, as our simulations were only set up to sim¬ 
ulate a single test particle. However, Figure [2] 
shows that there are several near-coplanar stable 
“islands” interior of Hydra, especially interior of 
Styx an d between Nix and K erberos. Observa¬ 
tions bv ISteffl fc SternI ( 2007h using HST placed 
very tight upper limits on the brightness of any 
rings in the Nix/Hydra (I/F<10“®), and therefore 
it is not likely that the region between Nix and 
Kerberos is populated by a dense ring. Thus the 
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most likely region for a ring to exist around Pluto 
is interior of Styx, just interior of the Charon 3:1 
MMR. 

5. Summary 

1. We determined new constraints for the 
masses of the known satellites of Pluto, in¬ 
cluding a lower mass for Nix, which allowed 
for Styx to have realistically-constrained 
mass. 

2. We found that additional low-inclination 
satellites are stable interior of Styx, between 
Styx and Nix, and between Nix and Ker¬ 
beros. Coorbitals of the known satellites are 
also possible, but must be very low mass. 

3. We found that retrograde orbits between the 
2:1 and 3:1 MMRs with Charon were quite 
stable, though very unlikely with current for¬ 
mation models of the Pluto-Charon system. 

This work was supported by the NASA New 
Horizons mission (contract number ^NAS02008). 
Special thanks to Marc Buie for providing the 
computing resources for this project, as well as 
Mark Showalter, Douglas Hamilton, and John 
Spencer for insightful discussions. 
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Fig. 1.— Derived mass uncertainties of the known 
Pluto satellites, uncertainty from observations and 
constrained to solutions stable for >10® Charon 
orbits. Histograms represent our normalized mass 
distribution for each satel lite, and the da s hed li nes 


indicate the masses from Brozovic et al. ( 2015f) . 



^ ^ ° Density of Orbits 

Average Barycentric Orbit Period (Charon Periods) 


Fig. 2.— Distribution of 10® possible low- 
inclination orbits. The blue points show all the 
orbits that were stable for >10® Charon periods. 
The red points are the 10% of the simulations with 
the lowest astrometric residuals, and so are the 
most likely orbits for an additional low-inclination 
satellite. The near-coplanar islands interior of 
Styx and between Nix and Kerberos are the most 

- stable regions for rings. 

This 2-column preprint was prepared with the AAS lAT^ 
macros v5.2. 
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Table 1 

Masses of Known Satellites (10“^ km^ s“^) 


Satellite 

Tholen et al. 
(2008) 

Youdin et al. 
(2012) 

Brozovic et al. 
(2015) 

This 

Work'^ 

Styx 

- 

- 

< 1.0 

0.10 ±0.03 

Nix 

39 ±34 

< 3.3 

3.0 ±2.7 

2.43 ±0.02 

Kerberos 

- 

- 

1.1 ±0.6 

1.02 ±0.05 

Hydra 

21 ±42 

< 6.0 

3.2 ±2.8 

3.69 ±0.04 


^1-17 uncertainty of solutions stable for >10^ Charon orbits 



i 5 o Relative Probability 

Average Barycentric Orbit Period (Charon Periods) 


Fig. 3.— Distribution of possible high-inclination 
orbits. The blue points show all the orbits that 
were stable for >10^ Charon periods. The red 
points are the 10% of the simulations with the 
lowest astrometric residuals. Retrograde-inclined 
orbits appear more stable simply because they in¬ 
teract less with the known satellites; such satellites 
could only plausibly be formed through a very low- 
probability capture event. 
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